Density functional study of Ni bulk, surfaces and the adsorbate systems Ni(lll) 
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Nickel bulk, the low index surfaces and the adsorbate systems Ni(lll) (\/3 x v / 3)R30°-Cl, and 
Ni(lll)(2 x 2)-K are studied with gradient corrected density functional calculations. It is demon- 
strated that an approach based on Gaussian type orbitals is capable of describing these systems. The 
preferred adsorption sites and geometries are in good agreement with the experiments. Compared 
to non-magnetic substrates, there does not appear to be a huge difference concerning the structural 
data and charge distribution. The magnetic moment of the nickel atoms closest to the adsorbate 
is reduced, and oscillations of the magnetic moments within the first few layers are observed in the 
case of chlorine as an adsorbate. The trends observed for the Mulliken populations of the adsorbates 
are consistent with changes in the core levels. 
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I. INTRODUCTION 



Adsorption on surfaces is one of the main focuses in 
surface science due to the enormous importance of topics 
such as corrosion or catalysis. A huge number of experi- 
ments has been performed to understand these processes. 
In addition, theoretical studies have become capable of 
q . modelling the adsorption process. 
Q ' The adsorption of alkali metals or haloges on metal- 
L ~ ~~ '. lie surfaces may be considered as two prototypes of ad- 
sorption processes. They are fairly simple and experi- 
mentally well studied. However, there have been quite a 
[* — ■ few surprises. For example, the occupancy of top sites 
by alkali metals as adsorbates was unexpected (for re- 
views, see references 1-3 and more general, reference 4). 
Halogens are another type of adsorbate which is exper- 
imentally well studied. On close-packed (111) surfaces, 
they usually occupy the face centered cubic (fee) hollow 
site, and in a few systems a partial occupancy (around 
15-25%) of hep (hexagonal close-packed) sites has also 
been observed 5,6 . In addition, there is now a growing 
I . set of simulations for halogens on metallic surfaces (e.g. 
chlorine 7-9 , bromine 10,11 or iodine 12 ). Both for alkali 
metals and for halogens as adsorbates, there are general 
questions such as the geometry, binding energies, diffu- 
sion barriers, charge transfer and the binding mechanism. 

Magnetic surfaces pose additional difficulties and ques- 
tions such as how the magnetic moment changes after 
adsorption. It is thus an interesting question which in- 
sights simulation can give for these systems. Nickel sur- 
faces and adsorption thereon have been well studied by 
experimental groups and thus are a prototype system 
for a computer simulation. For example, hydrogen on 
Ni(lll) 13 , benzene on Ni(lll) 14 , CH3 dehydrogenation 
on Ni(lll) 15 or CO on Ni(llO) 16 adsorption have been 
simulated recently. The present article will deal with the 
simulation of alkali and halogen adsorption, with CI and 
K as examples. 
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The intention of this article is to give a broad overview 
over nickel and the adsorbate systems. The computa- 
tional technique is based on a local basis set of Gaussian 
type functions. As nowadays most calculations, espe- 
cially for metals, rely on plane-waves, the following sec- 
tion II was devoted to the fairly simple systems of Ni bulk 
and the low index surfaces. This way, it will be demon- 
strated that results from a code based on Gaussian type 
orbitals reproduce the ones from plane- wave codes. In 
sections III and IV, the results for chlorine and potas- 
sium adsorption are presented, respectively. Finally, the 
results are summarized and discussed. Additional details 
on the calculations are given in the appendix. 

II. NI BULK AND LOW INDEX SURFACES 
A. Computational Parameters 

A local basis set formalism was employed where the ba- 
sis functions are Gaussian type orbitals centered at the 
atoms as implemented in the code CRYSTAL 17 . The cal- 
culations were done with spin-polarization, using the gra- 
dient corrected functional of Perdew, Burke and Ernzer- 
hof (PBE) 18 , and in some additional calculations the gra- 
dient corrected functional by Perdew and Wang (PW) 19 . 
These functionals perform actually nearly identical (see, 
e.g. references 18, 20, or with the CRYSTAL code for 
Ag 8 ), and also in the present article the results with both 
functionals are very similar. The PBE functional may be 
viewed as slightly favorable because it was designed to 
have a simpler form and derivation, and already revised 
versions exist (see, e.g., reference 20). 

The nickel basis set from reference 21 was employed, 
with the innermost [4s3pld] exponents being unchanged. 
In addition, 2 sp functions with exponent 0.63 and 0.13 
and a <i-function with exponent 0.38 as optimized in cal- 
culations for the bulk were used, so that the basis set 
as a whole was of the size [6s5p2d). The chlorine and 
potassium basis sets were the ones used in previous stud- 
ies on Cl/Cu(lll) and K/Cu(lll), respectively 7 ' 22 . For 
the calculations on the free atoms, enhanced basis sets 
with additional diffuse functions were used 23 . For the 
numerical integration, the potential was fitted with aux- 
iliary basis sets (for CI the same as in reference 7, for 
K as in reference 22, and for Ni the same as for Cu in 
reference 22). A fc-point net with 16 x 16 points in the 
surface Brillouin zone and 16 x 16 x 16 points for the 
bulk was used. A smearing temperature of 0.001 Eh was 
applied to the Fermi function (1^=27.2114 cV=315773 
K). This value for the temperature was chosen relatively 
low to ensure that the magnetic moment is not artificially 
modified by a too high value. For more details and tests 
on the computational parameters and basis sets, see also 
the appendix. 

The calculation of the charge and spin population is 
based on the Mulliken population, which is actually due 
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to McWeeny 24 ' 25 (see also reference 26). The Mullikcn 
charge is defined as P^S^^ (with density ma- 

trix P, overlap matrix S and A the set of basis functions 
for which the population is computed) . In the case of the 
spin population the density matrix for the corresponding 
spin polarization has to be used. This type of popula- 
tion analysis thus depends on the basis sets. It will be 
questionable if basis sets of different quality are used for 
various atoms: if one basis set is very small and the other 
very large, and especially if the larger basis set also has 
more diffuse exponents (i.e. smaller exponents, which 
results in a Gaussian with a relatively large spatial ex- 
tension), the larger basis set may describe a part of the 
charge of the atom with the smaller basis set as well, and 
thus the Mulliken charge will be artificially large for the 
atom with the larger basis set. It is thus important that 
balanced basis sets are used. Concerning the Mulliken 
populations computed in this article, the population can 
be expected to be fairly reliable in the case when iden- 
tical basis sets are used for the various atoms in a slab, 
so that the magnetic moments computed for the various 
layers of a nickel slab should not have a too large error 
due to the choice of the basis set. Also, comparing results 
for various adsorption sites should not be too unreason- 
able, because the basis sets for the systems are identical. 
This means, that, for example, a relative comparison of 
the Mulliken charge for chlorine on the fee site with the 
charge on the top site should not be too unrealistic, al- 
though the absolute value of the Mulliken charge may 
have a larger error. A more difficult issue is however, 
to compare the chlorine charge for Cl/Ni(lll) with the 
one for Cl/Ag(lll) or Cl/Cu(lll) because different basis 
sets for the metals are used and because of the different 
geometries. 

The adsorbate systems Ni(lll) (VSx V3)R30°-C1 and 
Ni(lll)(2 x 2)-K were modeled with a slab consisting of 
3 nickel layers and the adsorbate layer, with the adsor- 
bate layer on one side of the slab. The model is truly 
two-dimensional and not periodically repeated in the 
third dimension. This is different from the more com- 
mon approach where periodicity in three dimensions is 
used, and a number of vacuum layers is required to sep- 
arate the periodically repeated slabs. In the latter case, 
this also introduces an artificial electrical field due to 
the periodic boundary conditions in the third dimension, 
which is usually compensated for by an additional dipolc 
field 27 . The approach used by the CRYSTAL code for 
slabs avoids this artificial field as there are no periodic 
boundary conditions in the third dimension, and thus 
there is no correction required, and there is no parameter 
for the number of vacuum layers either. The simulation 
was performed for the fee, hep, bridge and top site (see 
figures 1 and 2). To ensure the stability of the results 
with respect to the number of layers, additional calcula- 
tions with thicker slabs were performed for the systems 
Cl/Ni(lll) (which is easier to converge than K/Ni(lll)). 
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B. Properties of bulk Ni and clean Ni surfaces 



Results for cohesive properties of bulk Ni are displayed 
in Table I. The computed values in this work repro- 
duce essentially the ones obtained with the gradient cor- 
rected functional of Perdew and Wang, with the plane- 
wave pseudopotential method 28 . The binding energy is 
overestimated, compared to the experiment. Also, the 
energy splitting between d s s 2 and cPs 1 is not well re- 
produced (the cPs 1 state should be lower by -0.03 eV, 
after averaging the experimental energies over spin-orbit 
components 29 ' 30 ). This is however, not a serious prob- 
lem, as the orbital occupancy in the solid is different 
from the atom and thus the error cancellation does not 
work too well, and a density functional description may 
be too crude to describe these effects. The splitting be- 
tween d 8 s 2 and dPs 1 state is a problem well known in 
quantum chemistry. The relativistic effect to the energy 
splitting between these states was computed to be of the 
order of 0.3 eV 30 ' 31 (i.e. in a non-relativistic calculation 
as the present one, the splitting should be -0.33 eV in- 
stead of -0.03 eV) and thus is not negligible. With highly 
accurate configuration interaction schemes 31 , the calcu- 
lated splitting is -0.18 eV; i.e. there is still an error of 
0.15 eV, compared to the experimental splitting, after 
relativistic effects were subtracted from the experimen- 
tal value. The value computed here (-0.035 E h = -0.95 
eV); is in reasonable agreement with other results from 
literature, e.g. in the range of -1.1 to -1.2 eV, both at the 
level of the local density approximation (LDA) and the 
level of the generalized gradient approximation (GGA), 
with the PW functional 28 . This difficulty in obtaining 
the correct energy splitting between d n ~ 1 s 1 and d n ~ 2 s 2 
was also identified as a general problem for the transition 
metals 32 , with the LDA, where a bias towards d n ~ 1 s 1 of 
~ 1 eV was found. With a different gradient corrected 
functional (B-LYP), an energy splitting of -0.82 eV was 
computed 33 . There is, however, no reason why other 
properties should be affected by this problem. 

The band structure is displayed in figure 3. It agrees 
well with the early tight binding 34,35 or augmented plane 
wave 36 results, and thus supports the validity of the ap- 
plied approach, relying on Gaussian type orbitals. 

Surface formation energies (in Eh per surface atom) 
of the clean, unrclaxed surfaces are displayed in Table 
II. The surfaces were modeled with a finite number of Ni 
layers (typically from 1 to 9). The PBE energies obtained 
here are within literature values. For the (111) surface, 
the surface energy was additionally computed with the 
Perdew- Wang functional, and is identical to the one with 
the PBE functional. 

Finally, in Table III, the magnetic moments of the lay- 
ers of the unrelaxed slabs are given. In all cases, they 
were computed using a 9 layer slab. They are in good 
agreement with literature values. For the (111) surface 
it is found that the second layer in the slab has the high- 
est magnetic moment and not the top layer, in agree- 
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ment with two of the published articles. For comparison, 
for the (111) surface, relaxation of the top layer was ad- 
ditionally taken into account. This resulted in a small 
inwards contraction by 0.03 A, which is consistent with 
experimental 37 and theoretical 38 findings. The magnetic 
moments were virtually unchanged, and thus relaxation 
does not seem to explain why in one publication it was 
observed that the outermost layer should have the highest 
moment 38 . For comparison, also the Perdew-Wang func- 
tional was used to compute the moments of the unrelaxed 
(111) surface, and the results of the PBE functional were 
essentially reproduced. As explained earlier, the calcula- 
tion of the magnetic moments is based on the Mullikcn 
population analysis. The error should, however, not be 
too large as identical basis sets are used for all the nickel 
atoms of the slabs. 



III. CL ON THE NI(lll) SURFACE 

As a first adsorbate system, the adsorption of CI on 
the Ni(lll) surface was investigated. This system was 
studied experimentally with LEED (low energy electron 
diffraction), Auger electron spectroscopy, work function 
measurements and flash desorption 39 , surface extended 
X-ray absorption fine structure (SEXAFS) and soft X-ray 
standing wave (XSW) 40 ' 41 ' 5 , and angle-resolved photoe- 
mission fine-structure measurements (ARPEFS) 42 . At 
low coverages there were indications of a (2 x 2) pat- 
tern, followed by a (V3 x \/3)R30 o structure, and a com- 
plex (11) pattern. The (V3 x V3)R30° structure is the 
one which has been extensively studied experimentally 
and thus also the calculations will focus on this pattern. 
There was agreement that the fee hollow was the pre- 
ferred site (with an fcc:hcp population ratio of around 
85:15 5 , and the Cl-Ni bond length was measured to be 
2.40 ± 0.04 A 40 , 2.33 ± 0.02 A 41 or 2.332 ± 0.006 A 42 . 
This corresponds to an interlayer distance of 1.92 A 40 , 
1.83 A 41 or 1.83 7 42 . The first nickel layer was found to 
be practically unrelaxed in two experiments 40,41 , whereas 
in the other article, a 5% contraction with respect to the 
bulk value was deduced 42 . 

It is thus an interesting issue to compare with results 
from the simulations. These results arc presented in Ta- 
ble IV. A structural optimization was performed with 
3 and 4 nickel layers. An additional single point cal- 
culation with 5 nickel layers was performed for fee and 
hep site, where the geometry from the thinner slabs was 
used, with one additional nickel layer at the distance cor- 
responding to the bulk. The vertical positions of the 
chlorine adatoms and of the first nickel layer were opti- 
mized. Lateral relaxations are not possible for fee and 
hep site. As bridge and top site are much higher in en- 
ergy, lateral relaxations were not taken into account for 
these sites either. Comparing the data for three and four 
layers, we see that the geometry is well converged, and 
the difference in the energy splitting between the sites 
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is within what is expected as the numerical noise (e.g. 
fee vs. hep site, the energy splitting is 0.7 mEh with 
three nickel layers, 0.2 mEh with four layers, and 0.6 
mEh with five layers). Firstly, we note that the fee hol- 
low is indeed found to be lowest in energy, with the hep 
hollow being nearly degenerate at an energy difference of 
only 0.2-0.7 mE h (~ 60-220 K). Although the numeri- 
cal noise is not negligible when energy differences on this 
small scale are computed (see appendix) , this energy dif- 
ference would be compatible with the simple thermody- 
namic estimate given in reference 6, with copper as sub- 
strate. The bridge site is higher in energy by ~ 3 mEh 
(or more than 900 K), the top site by 18 mEh, and thus 
thermodynamically fee and hep sites are clearly favored. 
The energy splitting between fee (or hep site), and the 
bridge site, may also be a crude estimate for the diffu- 
sion barrier from one threefold-hollow site to another one. 
These findings are actually very similar to those 7-9 for 
Cu(lll)(V3 x \/3)R30 o -Cl and Ag(lll)( v / 3 x V3)R30°- 
Cl, and also for the system Ag(lll)(l x l)-0 43 : on the 
(111) surface, the highly coordinated sites are clearly low- 
est in energy for adsorbates with a small radius. The 
bond lengths from chlorine to the nearest nickel dci-Ni nn 
also increase with increasing coordination number and 
the binding energy increases with increasing coordination 
number, as is to be expected from Pauling's argument 44 : 
one strong Cl-Ni bond in the case of the top site will 
be shorter than the three weaker bonds in the case of fee 
and hep, but the three bonds as a whole are energetically 
more favorable. The computed bond length agrees well 
with the values from the literature 40-42 . 

There are two different experimental results for the 
interlayer distance between first and second nickel layer. 
The computed value from this article agrees very well 
with the experimental from references 40, 41, whereas 
the large contraction of 0.10 A 42 with respect to the bulk 
appears to be unlikely. A smaller contraction would also 
be in line with simulations with Cu or Ag as substrate 7-9 . 

An effective radius for CI on Ni(lll) could be com- 
puted by subtracting an effective nickel radius of §4§A 

from the bond length rfci-Ni nn of 2.34 A. This would 
result in a radius of 1.09 A, in agreement with the com- 
puted values 7 ' 8 for CI on Cu(lll) (1.12 A) or CI on 
Ag(lll) (1.17 A). Thus the radius for CI is intermedi- 
ate between the atomic (0.99 A) and ionic (CI - ; 1.81 
A) 45 . The bond length d C \ -Ni nn is also in agreement 
with experimental values for systems with similar atomic 
radii for the substrate (e.g., bond lengths of 2.39 A for 
the substrate Cu(lll) 46 , 2.38 A for Rh(lll) 47 , 2.39 A 
for Pd(lll) 48 ). 

In Table V, the CI population is displayed. The data 
is extracted from the calculations where three nickel lay- 
ers were used to model the substrate, but it is virtually 
identical to the data obtained with four nickel layers. 
The negative charge is largest for the top site, smaller 
for the bridge site and smallest for hep and fee site. This 
is consistent with the position of the 3s core eigenvalue 
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relative to the Fermi level, which is lowest for hep and fee, 
increasing for bridge and largest for the top site (because 
the negative charge of the chlorine adsorbate is increas- 
ing and thus the core levels are destabilized). Note that 
the core level for majority spin is also slightly lower (by 
~ 0.002 Eh) than for minority spin. The density of states 
is visualized in figure 4, projected on the chlorine basis 
functions, and projected on all basis functions (i.e. the 
total density of states). The different occupancy of the 
two spin states becomes visible. Note that the chlorine 
spin is mainly due to the different population of the sec- 
ond and third p-shell, and not of the fourth (if we order 
the shells according to the size of the exponents, start- 
ing with the highest). This means that the outermost 
and thus most diffuse exponent does not carry a huge 
spin and thus the chlorine spin is not an artifact of the 
Mullikcn population (if the spin was on the most diffuse 
function only, one might argue that this basis function 
also describes the spin on the nickel atoms and thus the 
Mullikcn population analysis would be questionable). 

The computed Mullikcn charge is also in agreement 
with the experimental estimates of 0.1 to 0.2 |e| derived 
from work function measurements 39 . The overlap popu- 
lation X^eA X^es P^vSvn (with density matrix P, over- 
lap matrix S and A, B the set of basis functions of atoms 
A and B for which the overlap population is computed) 
is 0.09 for CI and the nearest nickel atoms. This may 
indicate that there is no evidence for a strong covalent 
contribution to the binding mechanism, similar to non- 
magnetic substrates 7 ' 8,22,49 , and now also in the case of 
a magnetic substrate. For comparison, the overlap popu- 
lation between two nearest nickel atoms in the top nickel 
layer is 0.10. However, this data should not be over- 
interpreted, and the overlap population should rather be 
viewed as a rough qualitative measure and not a quan- 
titative measure for covalency. Another indication that 
the covalency is not strong might be that in the case of 
the bridge site, where p x and p y orbitals do not have the 
same occupancy, the orbital with the larger overlap with 
the neighboring nickel atoms (p x in our choice of geome- 
try) has a slightly smaller population than p y (the total 
CI p population is 3.73 for p x , versus 3.82 for p y ). 

Charge density, charge density difference and spin den- 
sity are visualized in figure 5. As usual, the charge den- 
sity difference is more informative as the charge density. 
The transfer from the top nickel layer to the chlorine 
atom becomes obvious. In contrast to the charge density, 
the spin density also offers more information: it is very 
interesting that the spin density varies relatively strong 
still for the atoms in the third layer. This is confirmed 
by the data from Mulliken population analysis for the 
spin (table VI). The magnetic moments depend rela- 
tively strong on the number of layers, and the analysis 
has been performed by extracting the data from a single- 
point calculation for a slab with 5 nickel layers, with an 
identical geometry as for the 3-layer slab, i.e. the addi- 
tional layers were added at a distance as in the bulk. 

The data is presented in such a way that one number 
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is given, if all the atoms in a layer have the same distance 
to the chlorine atom, and two numbers, if there are two 
different distances for the atoms in one layer possible. 
For example, in the case of the fee site, looking at the 
atoms in the first or second layer, all three nickel atoms 
within a supercell have the same distance to the adsorbed 
chlorine. However, in the third layer there will be one 
atom which is closer than the other two. The population 
of this closer atom (0.59) is displayed in the left column 
of the results for the 3 rd layer, the data for the other two 
atoms (0.64) in the right column. 

It is interesting to note that those nickel atoms in the 
first layer which are closer to the chlorine atom have a 
reduced magnetic moment. This is strongest in the case 
of the top site where the atom vertically under the chlo- 
rine has its population reduced to 0.55, versus a moment 
of 0.76 for the other two atoms. It is also apparent for 
the bridge site (0.66 vs 0.73), and in the case of fee and 
hep site the data is identical because the nickel atoms in 
the first layer are identical by symmetry. 

In the second layer, the spin population is higher for 
those atoms which are closer to the adsorbate (this is 
possible for the hep and bridge structure). In the third 
layer, again those atoms which have a shorter distance 
to the adsorbate, have there moment reduced. These 
oscillations are still visible in the fourth and fifth layer. 

Although this data depends relatively strong on the 
number of layers, it becomes apparent that there are os- 
cillations in the spin density. It is also interesting to 
compare with the total charge: the latter quantity is vir- 
tually constant within the layers, from the second layer 
on, and has a value identical to that of a clean slab 
from the third layer on (also table VI). Only in the first 
layer, the atom(s) closer to the chlorine adsorbate have a 
slightly higher charge (i.e. more electronic charge) than 
the atom(s) which are further apart (visible for bridge 
and top site). 

Finally, the chlorine adsorbate has a magnetic moment 
of < 0.1 for all sites (the smallest value for the top site). 
The moment is always parallel to the nickel spin. 

It is interesting to note that in the simulations for CO 
on Ni(llO) 16 , it was also found that the moment of the 
atom(s) closest to the CO was reduced. 

IV. K ON THE NI(lll) SURFACE 

In addition, the system Ni(lll)(2 x 2)-K was studied 
with this approach. Experimentally, in a LEED study 
at various coverages, this was found to be the only com- 
mensurate structure 50 . In a further LEED study 51,52 , 
the top site was found to be occupied, with a K-Ni dis- 
tance of 2.82± 0.04 A, a vertical rumpling of the first 
Ni layer of 0.12 A, and a lateral displacement of 0.06 A 
of the three nickel atoms in the top layer which are not 
vertically under the K adsorbate. The adsorption site 
was later confirmed in SEXAFS measurements 53 , and a 
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larger bond length of 2.92± 0.02 A was obtained. With 
the same technique, also the bond length for the system 
K/Cu(lll) was computed and deduced to be by 0.06 ± 
0.04 A larger. From ARPEFS 54 , a top site adsorption 
with K-Ni bond-length of 3.02 ± 0.01 A was obtained. 
In addition, the first to second layer spacing was reduced 
by 0.13 A to 1.90 ± 0.04 A, but no substrate rumpling 
in the top layer was observed. Lateral displacements 
were not found either (0.00 ± 0.09 A) in this experiment. 
Finally, photoelectron diffraction (PhD) 55 resulted in a 
bond length of 2.87 ± 0.06 A, and a reduced first to sec- 
ond layer distance of 1.86 ± 0.06 A was observed, with 
virtually no rumpling in the first layer (0.01 ± 0.09 A). 

The experiments thus consistently support the top site 
as the adsorption site, but the geometrical parameters are 
not identical. Moreover, a comparison with the system 
K/Cu(lll) is very interesting, so that the influence of 
the partially occupied cZ-bands can be investigated. 

The results of the simulations are summarized in Ta- 
ble VII, with geometrical parameters as defined in figure 
6. Initially, the potassium position and the position of 
the top nickel layer was optimized (i.e. the uniform re- 
laxation of the top nickel layer). In a second step, a 
vertical rumpling of the nickel atoms in the top layer was 
allowed, as this is expected to be important for the top 
site occupancy. A further refinement would have been 
the additional possibility of lateral relaxations. These 
lateral relaxations were, however, found to be of minor 
importance for the related system K/Ag(lll), where the 
energy was lowered by the order of typically 0.1 - 0.2 
mEh, when these relaxations were allowed 49 . Thus, lat- 
eral relaxations were taken into account additionally only 
for the top site, which was found to be lowest in energy, 
and in order to compare with the experimental data for 
the lateral displacement. 

Firstly, the top site is found to be the most favor- 
able site. This is in agreement with the experimental 
findings 51-53 . By disabling the possibility of substrate 
rumpling, it is also demonstrated that substrate rumpling 
is crucial for top site occupancy, as it lowers the energy 
by 1.6 mEh, but only 0.4 mEh for the other sites. These 
findings are actually similar to the system K/Cu(lll) 22 . 
The top site thus becomes the most stable site because 
the nickel atom underneath is pushed into the substrate 
relative to its neighbors and the potassium with its rela- 
tively large radius thus also overlaps more with the sec- 
ond nearest nickel neighbors. 

The computed bond length (2.79 A) is also very sim- 
ilar to K/Cu(lll) (2.83 A) 22 so that these calculations 
do not indicate a major difference compared to Cu. Ex- 
perimentally, using the same technique for both adsor- 
bate systems, it was measured to be slightly shorter and 
it was speculated that this was due to greater adatom- 
substrate interactions involving the partially filled Ni d 
bands 53 . The possibility of a lateral displacement of the 
nickel atoms as observed experimentally (0.06 A) 51 was 
investigated. However, this could not be confirmed and 
the energy was found to be lowest for the case where this 
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lateral displacement is exactly zero. 

Bridge, fee and hep site are found to be energetically 
degenerate. The fact that all the sites are energetically 
very close (e.g. much closer than in the case of chlo- 
rine) is consistent with earlier findings 27,22,49 and was 
explained due to the large radius of the adsorbate which 
makes it experience only a small substrate electron den- 
sity corrugation 27 . This small energy splitting between 
the various sites is also consistent with the experiment 
because of various observations 53 : firstly, the SEXAFS 
data indicated enhanced in-plane motion resulting from 
a wide and shallow potential well. Secondly, no SEXAFS 
signal was observed for low coverages below 0.13 mono- 
layers. It was thus argued that lateral motion was es- 
sentially hindered by the presence of other adatoms. Fi- 
nally, with increasing temperature (already in the range 
of 120-145 K) the SEXAFS amplitude was found to de- 
crease. This was explained due to the occupancy of other 
sites, which leads to a destructive interference. 

The bond length c?K-Ni nn increases with the number 
of nearest nickel neighbors in the order top, bridge, fee 
and hep. The effective radius of the potassium adsorbate 
is 2.79-- 3 ^P A=1.54 A for the top site. This number is in 
good agreement with experimental values for potassium 
when occupying the top site 1 , or the computed value for 
K/Cu(lll) (1.55 A 22 ). 

Looking at the populations in Table V, we see that they 
(~ 0.3 |e|) are very similar for all considered sites, only for 
the top site the charge is slightly less positive. In agree- 
ment with this, the core eigenvalues of the 3s and 3p or- 
bital are virtually independent of the site. The core levels 
are identical for both spin polarizations which agrees with 
the finding that the magnetic moment of the potassium 
adatoms is negligible (see below). In the calculations for 
K/Ag(lll) 49 , it was recently shown that the adsorbate 
charge and thus also the core eigenvalue (and other prop- 
erties such as the bond length) change with the coverage. 
Such a change in the core eigenvalues was observed exper- 
imentally for alkali metals on metallic surfaces, e.g. for 
the systems Na/Cu(lll) and Na/Ni(lll) 56 , or for alkali 
metals on Ru(100) 57 or W(110) 58 . 

The overlap population is <~ 0.03 for K and the nearest 
nickel atom, so that there is also in this case no evidence 
for a strong covalent contribution to the binding mecha- 
nism. 

The density of states is displayed in figure 7. The inte- 
grated density of states, when projected on the potassium 
basis functions, differs so little for majority and minority 
spin that there is nearly no magnetic moment. For the 
top site, the spin is even antiparallel to the nickel sites, in 
contrast to the other sites (and in contrast to chlorine). 

Again, the charge density is not very informative, and 
it is more interesting to look at the charge density dif- 
ference (figure 8). Electrons have flown from the ad- 
sorbate mainly to the top nickel layer (and a relatively 
large part to the atom vertically under the potassium 
adatom). Thus, the charge of the potassium overlayer 
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has decreased. It is apparent that the spin density of 
the nickel atom vertically under the top atom is slightly 
lower than the spin density of its neighbors in the same 
layer, as the contours are closer to the nucleus. This is 
confirmed in the Mulliken spin population (table VIII). 
Again, those atoms in the top layer which are closest 
to the potassium adsorbate, have their moment reduced. 
No noteworthy variation is found in the second and third 
nickel layer, and the values for the spin population are 
very close to the ones for the clean Ni slab. The moment 
of the potassium adsorbate is negligible, as its magni- 
tude is < 0.003. Also, the charge only varies in the first 
layer, and already the second and third layer have charges 
which are identical to that of a clean slab. 

There are thus no oscillations of the spin population 
beyond the first nickel layer, probably because the potas- 
sium magnetic moment is much smaller than that of chlo- 
rine. 



V. SUMMARY 

We have studied Ni bulk, the clean low index sur- 
faces and the adsorption of chlorine and potassium on 
the Ni(lll) surface. 

The results for the adsorption geometry confirm the 
experimental findings. It turns out that there is no huge 
difference to copper as substrate; the preferred adsorp- 
tion site and energies are very similar. 

The reason for the preference of the top site for K is 
substrate rumpling which helps to increase the overlap of 
the K adsorbate with the second nearest nickel neighbors. 
This rumpling is thus crucial for the site preference. In 
the case of chlorine, the threefold hollow sites are very 
close in energy. In general, the adsorbate atoms prefer 
sites with the highest possible coordination number, or 
they try to increase the overlap to neighboring atoms by 
substrate rumpling. 

There appear to be no strong covalent contributions 
to the binding so that the mechanism appears to be 
mainly an ionic bond (although the charge transfer is 
rather small), and in the case of potassium, additionally 
a metallic binding mechanism. However, the quantitative 
analysis of the binding mechanism is certainly difficult. 
The computed charges seem to be fairly reliable, but the 
overlap population is probably only a qualitative tool to 
estimate the degree of covalency. 

Charges and core levels were found to be consistent, 
i.e. with increasing negative charge of the adsorbate, the 
core levels of the adsorbate are slightly destabilized. It is 
also demonstrated that the magnetic moments computed 
with the help of the Mulliken population analysis agree 
well with data from other schemes, e.g. pseudopotential 
plane wave. The calculation of work functions, however, 
seems to be difficult with a local basis set and the results 
are strongly basis-set dependent. 

Finally, the magnetic moment of those nickel atoms 
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closest to the adsorbate are always reduced. Oscillations 
of the spin population within one nickel layer are notable 
for the first three layers for Cl/Ni(lll), but only for the 
first nickel layer for K/Ni(lll), probably because of the 
very small magnetic moment of the potassium adsorbate, 
compared to chlorine. 
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APPENDIX A: TEST OF COMPUTATIONAL 
PARAMETERS 

The calculations require various sets of parameters, 
with the most important being: a finite number of k- 
points, a finite temperature to smoothen the integrand 
and thus to facilitate the numerical integration of the 
exchange-correlation functional, basis set parameters, 
and truncation schemes to reduce the number of matrix 
elements (i.e. integrals of operators such as kinetic en- 
ergy, nuclear attraction, coulomb and exchange, with the 
basis functions). In this first part of the appendix, the 
stability of the results with respect to the variation of 
parameters concerning the grid and the parameters con- 
cerning the selection of the integrals is tested. For further 
tests, see also the extensive tests performed for lithium 59 , 
or tests for copper 7 and silver 49 (fc-points, smearing tem- 
perature). In tables IX and X, the binding energies are 
computed with the default grid and truncation parame- 
ters, and with an improved grid and two sets of stricter 
and thus more accurate truncation parameters. 

The improved grid has roughly six times more sam- 
pling points than the default grid for the various struc- 
tures. As is displayed in the tables, the grid has virtually 
no impact on the energy splitting. For Cl/Ni(lll), the 
binding energy changes by at most 0.2 mEh for fee, hep 
and bridge site, and only for the top site an energy change 
of 0.9 mEh is observed, with respect to the default grid. 
This has, however, no impact as the top site is much 
higher in energy anyway. For K/Ni(lll) the absolute 
value of the binding energy is shifted by 0.7 mEh, and 
the relative energies of the various sites are identical. 

There are two truncation parameters in calculations 
when no Fock exchange is involved 17 : a general one 
("ITOL 1") which leads to a neglect of all integrals with 
an overlap below a certain threshold (default: 10~ 6 ), and 
a second parameter ("ITOL 2") which leads to the evalu- 
ation of the coulomb integrals at a lower level of accuracy 
by means of a multipolar expansion, if the overlap is be- 
low another threshold (default: 10 -6 ). As the results 
turned out to be more sensitive to these parameters, two 
sets of higher thresholds were used (10~ 7 for ITOL 1 and 



12 



2, and 1(T 8 for ITOL 1 and 2). In the case of chlo- 
rine, the largest error is found for the splitting between 
fee and hep site: it ranges from 0.7 to 2.1 mEh, for the 
various parameters, with the fee site always being more 
favorable. The bridge site is always clearly higher than 
the fee site by 3.3 to 4.2 mEh, and the top site is much 
higher in energy. For K/Ni(lll), the impact of these pa- 
rameters is smaller, and for example, the splitting from 
the top site to the site closest in energy ranges from 0.5 
mE h to 0.8 mE h . 

As a whole, the impact of these computational pa- 
rameters is certainly not negligible, but the conclusions 
are not affected. The fee site is favored for Cl/Ni(lll), 
with the hep site being nearly degenerate. In the case of 
K/Ni(lll), the top site is favored, with the other sites 
being very close in energy. 



APPENDIX B: BASIS SET DEPENDENCE OF 
THE RESULTS 

The choice of the basis set is of high importance for the 
quality of the results. In this part of the appendix, the 
results of various tests are displayed, to investigate the 
dependence of the results on the basis set. These tests 
focus on nickel bulk, the clean surfaces, and Cl/Ni(lll). 

Firstly, one might consider reoptimizing the exponents 
for the clean surface, for example keeping the exponents 
of the inner layer like in the bulk and reoptimizing the 
exponents of the outermost layers. This would result 
(with a slab with three layers) in exponents of 0.12 (sp, 
instead of 0.13 for the bulk) and 0.38 (d, as in bulk nickel) 
for the outer layers. This is thus a tiny change only and 
the results are virtually identical, for example the surface 
energy per atom changes by 0.6 mEh which is negligible. 
If an outermost sp-exponent of 0.12 instead of 0.13 is 
used for the bulk, exactly the same results as in table I 
are obtained. 

There is however a property which depends extremely 
strong on the basis set: the computed work function 
changes drastically when the outermost exponent is 
changed by a tiny value. For example, with the sp- 
exponent used for all the calculations (0.13), the work 
function of the Ni(100) surface is 0.142 E h (exp. 60 - 61 : 
0.192 E h ), of the Ni(110) surface 0.129 E h (exp. 60 ' 61 : 
0.185 E h ), and of the Ni(lll) surface 0.155 E h (exp. 60 - 61 : 
0.197 Eh). This data is practically stable with respect to 
the number of layers (i.e. 3 layers are sufficient to obtain 
a stable number). When, however, for example, the out- 
ermost sp-exponent of the outer layers is changed from 
0.13 to 0.12, the work function of the Ni(lll) surface 
changes by 0.01 Eh to 0.165 Eh, and it further increases 
with smaller exponents (e.g. 0.183 Eh, if the outermost 
sp-exponent of the outer layers is 0.10 instead of 0.13). 
Also, increasing this exponent leads to a lower work func- 
tion. The same problem shows also up with other metals, 
such as for example copper or silver (and also when the 
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exponents of all the layers and not just the outermost 
layers are changed). It appears thus that the work func- 
tion is a quantity which is not well described with a local 
basis set. The work functions may be qualitatively rea- 
sonable, but a quantitative comparison seems unadvised. 
Thus, also the work function data computed for Ag and 
the adsorption thereon 49 has to be considered with cau- 
tion. A problem with computing work functions with a 
local basis set was already discussed earlier 62 . 

The results thus indicate that the only property which 
appears to be strongly basis set dependent is the work 
function (total energy, surface energy, bulk modulus and 
lattice constant are essentially stable with respect to 
small variations in the basis set). Mulliken charge and 
magnetic moment are stable if an identical basis is used 
for all the nickel atoms in a slab. Of course, if basis sets 
of different type are used for the inner layer of a slab and 
the outermost layers, the Mulliken population also gets 
less reliable. 

Further tests were performed for the adsorbate sys- 
tem Cl/Ni(lll). One test was done by choosing a more 
tight outermost sp-exponent of 0.15 instead of 0.09 for 
chlorine. When a structural optimization is performed, 
virtually the same structural data as in table IV is ob- 
tained (the maximum deviation is 0.01 A). The binding 
energy is lower by 3mEh because of the poorer chlorine 
basis set. However, what is important is that the energy 
splitting for the various sites remains practically identical 
(the binding energy is -0.1326 Eh for the fee site, -0.1320 
Eh for the hep site, -0.1292 Eh for the bridge site, and 
-0.1151 Eh for the top site). Also, properties such as the 
population, magnetic moment and the relative position 
of the core levels with respect to the Fermi energy are 
essentially stable with respect to this change of the basis 
set (but again the work function changes strongly). 

As a whole, the basis set needs of course to be care- 
fully chosen. When this is done, the results for the ge- 
ometry, energies, core eigenvalues, populations, and mag- 
netic moments are reliable and not too strongly basis set 
dependent. The accurate determination of the work func- 
tion, however, appears to be difficult with a local basis 
set. 
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TABLE I. The ground state properties of bulk Ni. 





lattice constant do [A] E co h [Eh] 


B [GPa] magnetic moment [hb] 


PBE, this work 


3.53 0.218" 0.184 6 


203 


0.62 


PW, this work 


3.53 0.220" 0.185 6 


203 


0.61 


Ref. 28, LDA 


3.43 


255 


0.59 


Ref. 28, GGA 


3.53 


195 


0.61 


exp. 


3.52 63 0.163", 0.164 1 ' 63 


190 63 


0.61 45 




a energy with respect to a free nickel atom 


in its d s s 2 state 






b energy with respect to a free nickel atom 


in its cPs 1 state 






TABLE II. The surface energy [ sur f a <^ h atom ] of the low index nickel surfaces. 


surface 


PBE, this work Ref. 38 


Ref. 64 


Ref. 65 


(100) 


0.038 0.031 


0.039 


0.036 


(110) 


0.056 0.046 




0.049 


(111) 


0.028 (PW: 0.028) 0.024 


0.033 


0.026 


TABLE III. Magnetic moments [/ib] of the low index nickel surfaces, computed using a 


slab with 9 Ni layers. 


layer 


PBE, this work PW, this work Ref. 38 (9 layers, 






3 surface layers < 


m both Ref. 64 Ref. 66 Ref. 67 




sides are allowed to relax) 


(7 layers) (7 layers) 


(100) surface 








S 


0.729 0.76 


0.69 


0.68 


S-l 


0.634 0.68 


0.64 


0.60 


S-2 


0.632 0.66 


0.66 


0.59 


S-3 


0.611 


0.64 


0.56 


S-4 


0.619 






(110) surface 








S 


0.757 0.76 






S-l 


0.637 0.66 






S-2 


0.618 0.64 






S-3 


0.620 






S-4 


0.607 






(111) surface 








S 


0.648 (relaxed: 0.649) 0.638 0.68 


0.62 


0.63 


S-l 


0.658 (relaxed: 0.657) 0.649 0.65 


0.67 


0.64 


S-2 


0.622 (relaxed: 0.622) 0.613 0.62 


0.65 


0.58 


S-3 


0.617 (relaxed: 0.617) 0.604 


0.63 


0.58 



S-4 0.613 (relaxed: 0.613) 0.605 
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TABLE IV. Adsorption of CI on the Ni(lll) surface. c(ci-Ni top layer is the interlayer distance between the CI layer and the 
top Ni layer, dci-Ni nn is the bond length between CI and nearest neighbor Ni. djvii-jvi2 is the distance between first and second 
nickel layer. The distance between second and third nickel layer dNi2-Ni3 is held fixed at the bulk value (and for the thicker slabs 
also djvi3-jvi4 and d,Ni4-Ni5 are fixed at the bulk value). The adsorption energy is the difference E C i a t Ni(in) — ^Ni(in) — Eqi- 



Site 


dd— Ni top layer 


d,Nil-Ni2 


dci-Ni nn 


E adsorption 




[A] 


3 nickel layers 


[A] 


[A] 


r E h i 

L Cl atom J 


fee 


1.84 




2.02 


2.34 


-0.1358 


hep 


1.85 




2.02 


2.35 


-0.1351 


bridge 


1.90 




2.02 


2.27 


-0.1325 


top 


2.14 


4 nickel layers 


2.02 


2.14 


-0.1177 


fee 


1.83 




2.02 


2.33 


-0.1338 


hep 


1.84 




2.02 


2.34 


-0.1336 


bridge 


1.89 




2.02 


2.27 


-0.1307 


top 


2.14 


5 nickel layers 


2.02 


2.14 


-0.1159 


fee (single point calculation) 










-0.1337 


hep (single point calculation) 










-0.1331 



exp. (ARPEFS) 42 1.837 1.92 6 2.332 

exp. (SEXAFS, XSW) 40 ' 41 1.92, 1.83 2.03, 2.02 2.40, 2.33 



TABLE V. Charge and position of the 3s eigenvalue for Cl and the 3s and 3p eigenvalues for K on different adsorption sites, 
relative to the Fermi energy. For chlorine, the peaks are at slightly different positions for majority and minority bands, for 



potassium they 


are at the same 


position. 






site 


charge, in e 


Cl 


3s level, relative to Ef [E h ] 

on Ni: 


3p level, relative to Ef [Eh] 


fee 


-0.050 




-0.589 (majority spin); -0.586 (minority spin) 




hep 


-0.046 




-0.586 ; -0.584 




bridge 


-0.068 




-0.580 ; -0.578 




top 


-0.152 


K 


-0.537 ; -0.535 

on Ni: 




fee 


+0.294 




-1.178 (both spins) 


-0.583 (both spins) 


hep 


+0.295 




-1.178 


-0.582 


bridge 


+0.291 




-1.178 


-0.583 


top 


+0.269 




-1.179 


-0.583 
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TABLE VI. Mulliken spin and charge population for Cl/Ni(lll), on different adsorption sites, extracted from calculations 
with 5 nickel layers. The population for the nickel atoms in the individual layers is given. The number in parenthesis gives the 
number of atoms which are identical because of symmetry and thus have an identical population. The left column corresponds 
to that atom(s) in the layer which is closer to the chlorine adsorbate. For example, in the fee case, all three nickel atoms of a 
supercell in the first and second layer have the same distance and thus the population is identical, but one atom in the third 
layer is closer than the other two atoms resulting in two different values for the spin population. In the case of the clean slab, 
all atoms within one layer are degenerate and thus only one number is displayed. 

site I s * layer 2 nd layer 3 rd layer 4 th layer 5 th layer CI 

spin 

fee 0.67 (3 atoms) 0.64 (3 atoms) 0.59 (1 atom) 0.64 (2 atoms) 0.66 (3 atoms) 0.65 (3 atoms) 0.08 

hep 0.68 (3 atoms) 0.68 (1) 0.63 (2) 0.62 (3) 0.66 (3) 0.63 (1) 0.66 (2) 0.08 

bridge 0.66 (2) 0.73 (1) 0.67 (1) 0.63 (2) 0.60 (1) 0.63 (2) 0.66 (2) 0.65 (1) 0.64 (1) 0.66 (2) 0.08 

top 0.55 (1) 0.76 (2) 0.64 (3) 0.62 (3) 0.67 (1) 0.64 (2) 0.65 (3) 0.04 

clean slab 0.66 0.66 0.63 0.66 0.66 

charge 

fee 27.97 (3) 28.02 (3) 28.00 (1) 28.00 (2) 28.05 (3) 27.95 (3) 17.05 

hep 27.97 (3) 28.02 (1) 28.00 (2) 28.00 (3) 28.05 (3) 27.95 (1) 27.95 (2) 17.05 

bridge 28.00 (2) 27.90 (1) 28.02 (1) 28.01 (2) 28.00 (1) 28.00 (2) 28.05 (2) 28.05 (1) 27.95 (1) 27.95 (2) 17.07 

top 28.12 (1) 27.84 (2) 28.02 (3) 28.00 (3) 28.05 (1) 28.05 (2) 27.95 (3) 17.15 
clean slab 27.95 28.05 28.00 28.05 27.95 



TABLE VII. Adsorption of K on the Ni(lll) surface. cdt-Niia is the interlayer distance between the K layer and the layer 
made of those nickel atoms in the first layer closest to the K layer. diMiib-Ni2 is the distance between the layer made of those nickel 
atoms of the first layer which have moved away from the adsorbate and the second Ni layer. 8 = dNiia-Niib is the rumpling 
within the first layer. dn-Ni nn is the bond length between K and nearest neighbor Ni. Again, the distance between second and 
third nickel layer d Ni2 _Ni3 is held fixed at the bulk value. The adsorption energy is the difference E K at N i(iii) — ^Ni(in) — E K - 



Site 



d K - 



[A] 



d 



Nilb-Ni2 



[A] 



8 

[A] 



dn-i 



[A] 



*-> adsorption 

r E h i 

L K atom} 



fee 
hep 
bridge 
top 



2.71 
2.70 
2.70 
2.76 



without rumpling 

2.01 

2.01 

2.01 

2.01 



3.07 
3.06 
2.98 
2.76 



-0.0532 
-0.0531 
-0.0531 
-0.0528 



fee 
hep 
bridge 
top 



2.63 
2.63 
2.65 
2.68 



with rumpling 
1.99 
1.99 
1.97 
1.93 



+0.07 
+0.07 
+0.07 
+0.11 



3.06 
3.06 
2.99 
2.79 



-0.0536 
-0.0535 
-0.0535 
-0.0544 



exp. (LEED) 51 ' 52 
exp. (SEXAFS) 53 
exp. (ARPEFS) 54 
exp. (PhD) 55 



2.70+0.04 



experiment 
1.90+0.03 

1.90+0.04 
1.86 + 0.06 



0.12+0.02 

0.00+0.03 
0.01+0.09 



2.82+0.04 
2.92+0.02 
3.02 ± 0.01 
2.87+0.06 
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TABLE VIII. Mulliken spin and charge population for K/Ni(lll), on different adsorption sites, extracted from calculations 
with 3 nickel layers. The population for the nickel atoms in the individual layers is given. The left column corresponds to that 
atom in the layer which is closest to the potassium adsorbate. For example, in the top case, one atom in the first layer is closer 
than the other three atoms resulting in two different values for the population. In the case of the clean slab, all atoms within 
one layer are degenerate and thus only one number is displayed. Also, when the populations were found to be identical, only 
one number was displayed (usually from the second Ni layer on, marked with 'all'). 



site 



1 st layer 



2 nd layer 



3 rd layer 



K 



fee 0.58 (3 atoms) 0.62 (1 atom) 

hep 0.59 (3 atoms) 0.61 (1 atom) 

bridge 0.58 (2 atoms) 0.61 (2 atoms) 

top 0.55 (1 atom)0.61 (3 atoms) 
clean slab 0.66 



spin 

0.68 (all) 0.66 (all) 0.001 

0.68 (all) 0.66 (all) 0.001 

0.68 (all) 0.66 (all) 0.000 

0.68 (3 atoms) 0.69 (1 atom) 0.66 (all) -0.003 

0.69 0.66 



fee 28.04 (3 atoms)28.00 (1 atom) 

hep 28.04 (3 atoms)28.00 (1 atom) 

bridge 28.05 (2 atoms) 28.00 (2 atoms) 

top 28.08 (1 atom) 28.00 (3 atoms) 

clean slab 27.95 



charge 



28.10 (all) 
28.10 (all) 
28.10 (all) 
28.10 (all) 
28.09 



27.95 (all) 
27.95 (all) 
27.95(all) 
27.95 (all) 
27.95 



18.71 
18.71 
18.71 
18.73 



TABLE IX. Binding energies for Cl/Ni(lll), extracted from calculations with 3 nickel layers, as a function of the computa- 
tional parameters. The geometry was held fixed at the one optimized with the default parameters (table IV). 



site default grid, better grid, default grid, default grid, 

default ITOL (10~ 6 , 10~ 6 ) default ITOL (10~ 6 , 10~ 6 ) higher ITOL (10~ 7 , 10~ 7 ) even higher ITOL (10~ 8 , 10~ 8 ) 
fee -0.1358 -0.1357 -0.1359 -0.1357 
hep -0.1351 -0.1349 -0.1338 -0.1339 
bridge -0.1325 -0.1325 -0.1317 -0.1317 
top -0.1177 -0.1186 -0.1180 -0.1179 



TABLE X. Binding energies for K/Ni(lll), extracted from calculations with 3 nickel layers, as a function of the computational 
parameters. The geometry was held fixed at the one optimized with the default parameters (table VII). 



site default grid, better grid, default grid, default grid, 

default ITOL (10~ 6 , 10~ 6 ) default ITOL (10~ 6 , 10~ 6 ) higher ITOL (10~ 7 , 10~ 7 ) even higher ITOL (10~ 8 , 10~ 8 ) 
fee -0.0536 -0.0529 -0.0534 -0.0533 
hep -0.0535 -0.0528 -0.0532 -0.0531 
bridge -0.0535 -0.0528 -0.0532 -0.0532 
top -0.0544 -0.0537 -0.0539 -0.0541 
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FIG. 1. The threefold hollow structures considered for CI, adsorbed on the Ni(lll) surface, at a coverage of one third of a 
monolayer, (s/3 x y / 3)R30° unit cell. The nickel atoms in the top layer are displayed by open circles. The displayed chlorine 
adsorption sites are the threefold hollow sites (fee or hep hollow, filled circles). Note that these threefold hollow sites can not 
be distinguished in this figure. The ratio of the radii is chosen according to the computed values, i.e. using a nickel radius of 
1.25 A and a chlorine radius of 1.09 A. The top site would be vertically above a nickel atom, the bridge site vertically above a 
point which is in the middle between two neighboring nickel atoms. 




X 
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FIG. 2. Potassium, adsorbed on the Ni(lll) surface, at a coverage of one fourth of a monolayer, (2 x 2) unit cell. The nickel 
atoms in the top layer are displayed by open circles. The displayed potassium adsorption site is the top site vertically above 
the nickel atoms (filled circles). The ratio of the radii is chosen according to the computed values, i.e. using a nickel radius of 
1.25 A and a potassium radius of 1.54 A. 
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FIG. 3. Band structure of ferromagnetic nickel, with the PBE functional. The majority bands are drawn with full lines, the 
minority bands with dashed lines. 
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FIG. 4. The density of states, projected on the chlorine basis functions, (upper two graphs); and projected on all basis 
functions (lower two graphs), for majority (first and third graph) and minority band (second and fourth graph). 



25 



FIG. 5. Charge density (left), charge density difference (middle) and spin density (right) for CI on Ni(lll), from a three-layer 
slab, fee site. The charge density is displayed with full lines, increasing in steps of 0.01, from onwards. The charge density 
difference (charge of adsorbate system minus charge of a layer of chlorine atoms minus charge of a clean nickel surface) is in 
steps of 0.001, excess negative charge is displayed with full lines (from onwards), electron depletion with dashed lines (from 
-0.001 onwards). Positive spin density is displayed with full lines, increasing from in steps of 0.001. Negative spin density is 
displayed with dashed lines, decreasing from -0.0002 in steps of -0.0002. The plane goes through the center of the CI adsorbate, 
and through the nickel atom in the third layer vertically under the chlorine, and through the nearest neighbors of this nickel 
atom in the third layer. 
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FIG. 6. Geometrical parameters for the adsorption studies. 
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FIG. 7. The density of states, projected on the potassium basis functions, (upper two graphs); and projected on all basis 
functions (lower two graphs), for majority (first and third graph) and minority band (second and fourth graph). 
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FIG. 8. Charge density (left), charge density difference (middle) and spin density (right) for K on Ni(lll), from a three-layer 
slab, top site. The charge density is displayed with full lines, increasing in steps of 0.01, from onwards. The charge density 
difference (charge of adsorbate system minus charge of a layer of potassium atoms minus charge of a clean nickel surface) is in 
steps of 0.001, excess negative charge is displayed with full lines (from onwards), electron depletion with dashed lines (from 
-0.001 onwards). Positive spin density is displayed with full lines, increasing from in steps of 0.001. Negative spin density is 
displayed with dashed lines, decreasing from -0.0002 in steps of -0.0002. The plane goes through the center of the K adsorbate, 
and through the nickel atom vertically below, and through the nearest neighbors of this nickel atom in the first layer. 
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